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SUMMARY 


A  study  was  undertaken  to  examine  present  and  proposed  methods  of 
calculating  the  attenuation  afforded  "in  and  down"  scattered  radiation  by  a  horizontal 
slab,  and  the  effects  of  finite  fields  of  contamination  in  order  to  recommend  improve¬ 
ments  and  up  dating  for  OCD  documents  TR-20,  Volumes  I  and  2,  "Shelter  Design 
and  Analysis." 

In  this  study  we  examined  the  assumption  used  in  computing  the  "in  and 
down"  attenuation  factor,  assessed  their  effect  on  the  resulting  attenuation  factors, 
and  compared  the  results  with  the  latest  available  experimental  data.  Recommendations 
are  provided  as  to  "best  value"  factors  and  how  they  might  be  included  in  existing 
publications. 

The  effects  of  finite  fields  of  contamination  were  subjected  to  further 
analysis  by  directly  comparing  the  existing  method  of  calculation  with  results 
obtained  using  transmission  coefficients  calculated  by  the  Monte  Carlo  method.  This 
comparison,  though  complete  only  for  above-ground  locations,  indicates  that  further 
study  is  required  before  the  existing  method  of  calculation  can  be  relied  upon. 

The  study  of  in  and  down  scattered  radiation  is  described  in  Chapters  1 
through  3:  finite  fields  of  contamination  are  in  Chapter  4.  These  major  conclusions 
may  be  drcwn  from  this  study: 

The  "In  and  Down"  Attenuation  Factor 


1.  The  in  and  down  attenuation  factor,  (B'0(X,  u),  represents  a  significantly 
better  estimate  of  this  effect  than  the  values  previously  used  (Bg  (X)  ). 

2.  Calculated  and  measured  attenuation  factors  agree  with:n  18.5  per  cent; 
the  analytical  value  is  conservative.  It  was  found  that  the  standard 
deviation  of  the  ratio  of  the  attenuation  calculated,  Bg(X,u),  to  that 
measured  in  346  experimental  measurements,  was  13.5  per  cent. 


Finite  Fields  of  Contamination 


3.  The  current  method  of  computing  the  effects  of  finite  fields  of 
contamination  requires  further  study. 

20 

In  total  dose,  the  NBS  Monte  Carlo  compilation  method  and 
that  presently  proposed  ^  agree  fairly  well  for  detector  locations 
on  the  first  floor  of  a  sheltering  structure. 
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5.  Major  differences  in  the  fraction  of  scatter  to  direct  dose  were 
found  when  using  the  NBS  method,  casting  serious  doubt  on  the 
accuracy  of  current  methods  of  shelter  analysis  when  the  detector 
is  shielded  from  direct  radiation  (i.e.,  in  a  basementor  an 
upper-story  location  shadowed  by  a  heavy  floor). 
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CHAPTER  I 

IN  AND  DOWN  SCATTERED  RADIATION  -  THEORY 


INTRODUCTION 

Most  of  the  error  observed  between  measured  and  computed  values  of 
dose  :ote  »n  basement  areas  has  been  attributed  to  the  attenuation  provided  by  the 
basement  ceiling.  Values  of  this  attenuation  previously  were  computed  in  approxi¬ 
mate  form  for  cobalt,  cerium,  and  the  1.12  hr.  fallout  energy  spectrum,  using  a 
single  Compton  scatter  energy  change  and  transmission  coefficients  computed  by 
the  Monte  Carlo  method,  assuming  that  all  incident  radiation  lay  in  the  plane 
mutually  perpendicular  to  the  exterior  wall  and  the  ground  and  that  it  passed  through 
the  detector  location.  Agreement  between  the  attenuation  computed  in  this 
approximate  fashion  and  that  measured  was  much  better  than  that  obtained  with  the 
engineering  method. 

Our  purpose  in  this  study  has  been  to  re-examine  the  effects  of  the 
assumptions  used  and  to  recalculate  the  basement  ceiling  attenuation  more  accurately 
by  (1)  relaxing  the  restriction  that  incident  radiation  must  lie  in  the  plane  containing 
the  detector  that  is  mutually  perpendicular  to  the  ground  and  the  wall  and  (2)  using 
a  finer  gradient  when  performing  the  numerical  integration  over  the  contaminated  area. 


Basement  Ceiling  Attenuation 

Approximately  80  per  cent  of  the  radiation  arriving  at  the  exterior  wall 
of  a  structure  located  in  an  infinite  contaminated  field  is  uncollided.  The  remaining 
20  per  cent  is  singly  or  multiple  scattered  radiation  of  lower  energy.  This  "atmospheric" 
scattered  radiation  is  divided  about  evenly  between  the  upper  and  lower  hemisphere 
for  a  detector  located  above  ground  level.  Thus,  in  estimating  the  attenuation  afforded 
by  a  basement  ceilings  the  first  assumption  is  that  all  radiation  arriving  at  the 
above-ground  v-alis  of  a  structure  that  scatters  to  the  basement  area  is  uncollided. 
According  to  this  assumption,  all  scattered  radiation  is  of  higher  energy,  therefore 
rttenuating  less  rapidly  tnan  actual  scattered  radiation,  hence  producing  a  conservative 
estimate  of  the  radiation  scattered  to  a  below-ground  area.  To  simplify  the  problem, 
we  adopted  a  second  assumption:  radiation  scattered  by  the  above-ground  portion 
of  the  structure  is  reduced  in  energy  by  a  single  Compton  scatter,  through  the  angle 
required  to  reach  the  detector  only.  This  assumption  is  also  conservative  in  that  it 
treats  all  multiply  scattered  photons  as  of  higher  than  their  actual  energy.  Raso  has 
computed  the  average  energy  of  scattered  photons  for  incident  energy  of  1.25  mev 
and  finds  that  this  figure  is  about  85  per  cent  of  that  predicted  by  single-scotter 
theory.  The  third  major  assumption  is  that  the  detector  is  relatively  near  the  basement 


t 
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ceiling;  i.e.,  radiation  arriving  at  a  basement  detector  by  a  ceiling  slab  is  attenuated 
exactly  as  much  as  that  of  parallel  monoenergetic  radiation  of  the  same  energy  and 
angular  incidence  in  relation  to  the  external  surface  of  the  basement  ceiling. 

The  procedure  used  to  compute  the  ceiling  attenuation  was  as  follows 
(see  Figure  1).  The  contaminated  field  was  subdivided  into  "n"  increments  of  equal 
radiation  worth,  for  each  of  the  nine  energy  components  of  the  1.12  hr.  fallout 
spectrum,  when  viewed  by  a  detector  located  at  the  scattering  area.  The  total  angle 
of  scatter  (\jj)  then  was  determined  at  the  scattering  location  between  the  source 
patch  and  the  detector.  Notice  that  the  previously  used  model  for  computation 
was  a  cylinder  (or  the  equivalent,  the  short  wall  of  a  building  with  large  eccentricity) 
and  therefore  the  exit  azimuthal  angle  (3  was  identically  zero.  Then,  the  new 
scattered  energy  was  computed  from  the  Compton  scatter  equations.  The  dose 
arriving  at  the  detector  location  is  taken  as  proportional  to  the  amount  of  radiation 
energy  arriving  at  the  wall  from  the  source  patch,  reduced  by  the  ratio  of  the  initial 
energy  to  that  found  after  scattering,  and  multiplied  by  the  attenuation  introduced 
by  the  basement  ceiling  (taken  as  equal  to  that  of  parallel  monoenergetic  radiation 
of  identical  energy  and  angular  incidence)  and  the  air  absorption  coefficient. 

A  quantity  proportional  to  the  total  overall  dose  may  then  be  calculated  by  summing 
this  quantity  over  all  energies  of  the  fallout  spectrum  weighed  by  the  energy  spectrum 
and  over  the  infinite  plane  contaminated  area.  The  ceiling  attenuation  factor  was 
then  computed  by  dividing  this  quantity  by  the  similar  quantity  for  zero  ceiling 
thickness. 


This  may  be  expressed  in  equations  in  the  form  presented  here  (see 
Figure  1).  The  contaminated  field  beyond  the  wall  is  subdivided  into  a  total  of  "n" 
contaminated  patches  of  equal  worth  (when  viewed  by  a  detector  located  at  the  wall 
scattering  volume)  in  cylindrical  geometry.  The  total  uncollided  dose  rate  at  height  h 
from  a  semi -infinite  field  of  contamination  is: 


D(E)  =  const 


.-|it(E)p2,pdp 


const  Ej(p (E)h) 


.where: 


Mt(E) 


P 

E, 


=  is  the  total  coefficient  in  air  for  photon  energy  E 
=  the  height  of  the  scattering  element 
=  the  slant  distance  from  detector  to  differential  source  area 
=  the  exponent  integral  of  the  first  kind 
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The  total  number  of  differential  patches  "n"  is  composed  of  "i"  rings  of 
equal  worth,  each  subdivided  into  "j"  azimuthal  patches.  The  center  of  each  patch 
thus  is  specified  by  coordinates  9|  and  a^  relative  to  the  scattering  area. 


a 

m 


=  the  dzimuthal  coordinate  (see  Figure  1)  of  source  patch  "m" 


a 

m 


it(m  +  1/2) 
1 


the  polar  coordinate  of  a  source  patch  for  contamination  of 
ehergy  E 


H  t(E)h 

0j(E)  =  cos"1  Ht(Srp^ 


E^t(E)P;)  = 


E^M  f(E)h) 


T 

I 


tt+  1/2) 


For  heights  of  the  scattering  volume  "h"  that  are  small  with  respect  to 
a  mean  free  path,  the  quantity  9.(E)  is  nearly  independent  of  "h". 

I  ' 

The  detector  is  located  with  respect  to  the  scattering  volume  by  azimuthal 
angle  (3  and  polar  angle  9.  The  cosine  of  the  total  angle  through  which  a  photon  must 
scatter  from  the  source  patch  "I",  "m",  to  reach  the  detector  therefore  is: 


cos  4-' 


scos  9  cos  9.  -  sin  9  sin  9ft  cos  (180  -  a 

A  {  X  m 


which  is  implicitly  a  function  of  the  initial  energy  E0.  Thus  the  energy  E  in  mev  of 
the  radiation  scattered  b>  the  wall  toward  the  detector,  after  suffering  a  single  Compton 
scatter,  would  be 


E  =  E 

o 


.TRT 


.51 

TT 


cos  40 


; ' 

•j 
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If  the  probability  that  a  photon  will  scatter  from  source  patch  9|,a^  to 
a  detector  located  at  9,  (3  relative  to  the  scattering  volume,  is  written  as 
P(0|,  Qm  -'*■9,(3),  the  dose  seen  by  such  a  detector  located  below  a  floor  of  thickness  X 
from  a  single  Compton  scatter  is 


D(Eq,X,9,(3) 


A£  p{9„  a  9,  (3)  E  pn  (E)  cos  9  S'(X,  cos  9,E) 
o  l  m  u 

n  E 

o 


where: 


AE  =  the  total  quantity  of  energy  Eq  incident  on  the  scattering 
volume 

Eq  =  the  energy  of  interest 

E  =  the  energy  of  the  photon  after  suffering  a  single  scatter  to  the 

detector 

p^(E)  =  tissue  equivalent  energy  absorption  cross  section 

cos  9  S'(X,  cos  9,  uE)  =  attenuation  factor  for  parallel  monoenergetic 

radiation  of  energy  E,  angle  of  incidence  9, 
and  slab  thickness  X 

The  attenuation  introduced  by  the  ceiling  above  the  detector  therefore 
is  equal  to  the  integral  of  this  quantity  over  the  total  scattering  volume  (the  wall) 
divided  by  the  same  quantity,  with  thickness  X  equal  to  zero.  In  practice  the 
quantity  P(9|,  _»9, (3)  is  nearly  constant  and  may  be  removed  from  under  the 

integral  sign  and  canceled.  If  P(C|,  c^-^9, (3)  is  taken  as  approximately  equal  '■o 
the  differential  Compton  scatter  term  rather  than  constant,  the  differences  are  quite 
small.  This  is  Plustrated  in  Table  1  for  the  region  of  most  interest,  cos  9.  Here 
agreement  is  found  to  be  well  within  the  region  of  other  uncertainties  of  the  calculation. 
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TABLE  1 


Comparison  of  attenuation  factors  computed  with  probability  of  scatter 

equal  to  the  Compton  differential  scatter  probability 

and  a  constant  for  (3  *  0  and  E  *  1.25  mev 

o 


cos  0 

1.0 

.75 

.50 

.25 

_ _ _ 

.10 

Scatter 

probability 

Comp. 

Const. 

jjjjjjjj 

Const. 

mu 

Xf 

10  psf 

.90 

.89 

.83 

.82 

.68 

.67. 

.39 

.36 

.13 

.12 

20  psf 

.76 

.72 

.64 

.62 

.45 

.42 

.18 

.16 

.035 

.032 

40  psf 

.49 

.43 

.37 

.38 

.19 

.18 

.052 

.050 

.010 

.0091 

80  psf 

.18 

.12 

.11 

.029 

.039 

.038 

.0086 

.0070 

.0018 

Bb 

160  psf 

.023 

.010 

.0098 

.0068 

.018 

.0014 

.0038 

.00030 

.00081 

ror  simplicity,  the  scattering  probability  may  therefore  be  taken  as  constant 
over  the  range  of  interest  of  this  problem,  and  the  attenuation  at  position  0,  may 
be  written  as 


Att  (9,  h,  X,  Eo) 


\  p  (E)  cos  OS'fo  cos  0,E) 
c  a 
o 


|  H  (E)  cos  9S'(X  =  0,  cos  0, E) 
c  a 
o 


Two  aspects  of  the  problem  may  now  be  evaluated  quite  easily.,.  First, 
since  Raso^  found  that  the  average  energy  of  scattered  photons  was  about  85  per  cent 
of  Compton  single-scatter  energy,  what  effect  would  this  degradation  of  energy  of 
the  scattered  photons  have?  Second,  what  effect  other  than  the  straightforward 
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geometric  difference  of  slant  penetration  distance  in  the  floor  would  a  variation 
in  angle  beta  have  on  the  attenuation  factor? 

These  effects  are  illustrated  in  Table  2.  Here  again  it  is  evident  that 
essentially  there  is  no  secondary  geometric  effect  on  the  floor  attenuation  and  that 
the  assumption  of  scattered  energy  equal  to  85  per  cent  of  Compton  single  scatter 
lowers  the  attenuation  factors  by  only  a  small  amount. 

Since  most  buildings  are  other  than  cylindrical  in  shape  the  necessary 
integration  over  the  variable  beta  must  be  examined.  Previously,  it  was  assumed 
that  the  angular  distribution  of  radiation  emerging  from  the  rear  of  a  vertical  wall 
in  the  plane  parallel  te>  the  ground  is  closely  approximated  by  the  cosine  of  the 
emergent  angle  relative  to  a  perpendicular  to  the  scattering  surface.  It  was  demonstrated 
above  (see  Table  2)  that  the  principal  effect  of  this  angle  being  other  than  zero  is 
in  the  straightforward  geometric  increase  in  the  slant  penetration  distance  of  the 
slab  over  the  detector  that  a  photon  must  traverse  to  reach  a  detector  a  given  distance 
from  the  wall.  This  distance,  represented  equivalently  by  an  increase  in  the  angle  9 
to  91  for  an  attenuating  slab  of  constant  thickness,  may  be  calculated  quickly  by  the 
relationshi  p: 


Tan  9‘  =  tan  9/cos  |3 

The  attenuation  afforded  by  a  rectangular  structure  may,  then,  be  calculated  in  a 
manner  similar  to  that  applied  to  a  cylindrical  structure  by  integrating  the  dose 
scattered  from  a  given  scattering  volume  (represented  by  angle  |3),  weighed  by  the 
emergent  angular  distribution  for  a  finite  thickness  slab  divided  by  rhe  similar  quantity 
for  a  slab  of  zero  thickness.  Thus 


- 


all  walls 

i 

i 


AD  (91,  |3,X,  E  )  cos  (3d  (3 

o 


all  walls 


AD  (9’,  (3,  X  =  0,  E  )  cos  (3d  (3 
o 


where  8^(X,  Q,  Ep,  ^  )  -  the  attenuation  that  would  be  afforded  by  a  slab  to  in  and 
down  scattered  radiation  in  a  rectangular  structure. 
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L/W  =  length -to- width  ratio  of  the  structure 

This  quantity  has  been  calculated  for  several  conditions  of  floor  mass, 
building  size,  and  eccentricity  by  numerically  summing  over  ^he  azimuthal  dis¬ 
tribution  in  5°  increments.  The  results  of  this  series  of  computation  ure  shown 
in  Table  3,  where  they  are  compared  with  the  values  obtained  by  computing  the 
floor  attenuation  of  a  cylindrical  structure  of  equal  solid  angle  fraction  when 

TABLE  3 

Comparison  of  attenuation  factors  computed  in  circular 
and  rectangular  geometry  for  equal  sol  id  angle  fraction,  u 
E  =  1.25  mev 

(Note  Structure  is  of  differential  wall  height) 


viewed  by  the  detector.  This  table  indicates  that  these  two  methods  give  excellent 
agreement  and  that  there  is  therefore  no  significant  "geometry  factor"  for  floor 
attenuation. 

The  assumptions  made  in  calculating  the  attenuation  afforded  to  in  and 
down  scattered  radiation  by  a  floor  slab  upon  evaluation,  have  been  found  to 
effect  but  Sidle  the  predicted  attenuation  values.  Using  a  constant-scatter  probabili 
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rather  than  the  Compton  single-scatter  probability  slightly  decreases  the  predicted 
basement  detector  dose  rates,  whereas  assuming  scatter  energy  equal  to  the  single 
Compton  scatter  energy  and  neglecting  secondary  geometric  effects  in  the  rectangular 
wall  case  slightly  increases  predicted  dose  rate. 

The  floor  barrier  computation  as  described  above  has  been  performed 
for  the  approximate  energy  of  cesium  and  cobalt  radiation  as  well  as  for  each  of 
the  nine  energy  components  of  the  1.12  hr.  fallout  spectrum.  Figures  2  through  4 
present  the  results  of  this  series  of  computations  as  a  function  of  the  mass  thickness 
of  the  attenuating  slab  and  the  solid  angle  to  the  scattering  volume,  in  Figure  4 
is  presented  the  result  of  summing  the  data  for  the  nine  energy  groups  of  the  1.12  hr. 
fallout  spectrum  weighed  by  the  value  of  each  of  these  groups.  Previously  estimated 
values  produced  by  a  rough  hand  calculation  are  generally  within  about  10  per  cent 
of  the  values  shown. 

The  use  of  these  figures  is  quite  straightforward.  In  calculating  the  exact 
dose  in  the  basement  area,  the  dose  from  each  differential  increment  of  height  of 
the  scattering  wall  should  be  multiplied  by  the  attenuation  factor  (Figures  2  through 
4)  for  the  basement  ceiling  at  the  solid  angle  fraction  of  that  differential  wall  height 
(when  viewed  by  the  detector)  and  these  results  summed  over  the  entire  scattering 
wall  height.  It  was  demonstrated  earlier,  however,  that  an  equivalent  value, 
generally  within  10  per  cent  of  that  calculated  by  this  method,  may  be  obtained 
easily  by  averaging  the  solid  angle  to  the  top  and  bottom  of  the  scattering  wall  and 
determining  the  attenuation  factor  for  this  average  solid  angle  fraction.  Several 
examples  of  this  procedure  (for  a  structure  of  10  ft.  wall  height  with  a  detector  7  ft. 
below  the  floor)  are: 

h 


Bui  Iding  plan  area 

Wall  thickness 

Bq(u,X)* 

r  b  (u,x)dh 

J  o 

0 

20  x  20 

25 

0.45 

0.48 

50 

0.19 

0.20 

o 

X 

o 

25 

0.26 

0.28 

50 

0.086 

0.098 

100  X  100 

25 

0.080 

0.083 

50 

0.018 

0.019 

*More  properly,  these  quantities  should  be  B0(u,X)L(d)  and  /  B  (i*>,X)l(d)dh,  where 

"o  ° 

L(d)  represents  the  dose  variation  with  height  d  above  the  contaminated  plane.  If 
this  further  refinement  is  included,  both  quantities  are  reduced  slightly,  the  integral 
quantity  always  being  reduced  more  than  the  multiplication  quantity.  If  the  scattering 
wall  extends  from  zero  to  10  ft.,  the  integral  quantity  is  reduced  by  approximately 
8  per  cent,  while  B0(w,  X)L(d)  remains  essentially  the  same. 
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CHAPTER  2 

IN  AND  DOWN  SCATTERED  RADIATION  -  EXPERIMENT 

INTRODUCTION 


Over  the  past  few  years  several  experiments  have  provided  data  useful 
in  verifying  the  basement  ceiling  barrier  factor,  Bo(Xf,u  ),  developed  in  this  report^ 
Recent  experiments  at  the  Defense  Chemical,  Biological  and  Radiation  Laboratories 
in  Canada,  at  the  Protective  Structures  Development  Center, *“*  and  at  Technical 
Operations  Research,  '  were  performed  in  such  a  way  as  to  determine  Bg(Xf,u) 
by  comparing  measured  results  with  and  without  a  basement  ceiling  slab  in  position. 
In  addition,  several  other  total  building  experiments  have  been  performed  that  are 
useful  in  comparing  measured  and  computed  basement  data.  The  latter  experiments 
either  were  not  specifically  designed  for  "in  and  down"  measurements  or  the  only 
results  available  to  date  are  cf  a  preliminary  nature  limiting  their  comparison  to 
the  total  basement  protection  factor,  without  isolating  the  ceiling  barrier  factor 
component,  Bo(Xf,u).  Brief  comparisons  are  here  presented,  however,  between 
buiiding  attenuation  values  computed  by  means  of  the  new  barrier  curves  B^Xpu) 
of  this  report  and  the  measured  results.  Of  the  work  cited,  only  the  DCBRL  experi¬ 
ment  is  covered  in  detail  in  this  report.  The  remaining  experiments  have  been 
covered  in  the  Conesco  report  entitled  "The  Preparation  of  Simplified  Manuals 
for  Shielding  Analysis."^ 


"In  and  Down"  Scattering  Measutcments  at  the  Defense  Chemical,  Biological  and 
Radiation  Laboratories11  5 

The  experimental  arrangement  (see  Figure  5)  for  the  "In  and  Down" 
scattering  study  at  DCBRL^  consisted  of  measuring  the  dose  distributions  produced 
by  the  radiation  transmitted  and  scattered  from  a  vertical  barrier  that  penetrated  an 
adjacent  horizontal  barrier.  A  Cesium-137  source  collimated  to  a  10°  half  angle  by 
a  lead-filled  collimator  was  utilized  as  the  radiation  source.  This  radiation  beam 
impinged  on  a  3  ft.  high  by  4  ft.  long  vertical  concrete  slab,  providing  an  effective 
circular  pattern  of  incident  dose  spread  over  an  area  20.5  in.  in  radius.  This  area 
was  centered  at  the  physical  center  of  the  slab  on  the  slab  face.  Measurements  of 
the  scattered  dose  were  taken  along  a  vertical  traverse  positioned  12,  24,  and  36  inches 
to  the  rear  of  the  vertical  slab  at  depths  below  the  top  surface  of  the  horizontal  slab 
ranging  from  3-3/4  to  30  in.  Veit'cal  slab  thickness  variations  were  0  p>f,  6.45  psf 
plywood,  and  34  and  68  psf  concrete.  The  horizontal  slabs  were  of  0,  36.5,  and  73  psf 
r>  reinforced  concrete.  The  detectors  were  partially  shielded  with  lead  to  limit 
their  exposure  to  multiple  scattering  from  the  ground. 
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The  attenuation  attributed  to  the  horizontal  (ceiling)  barrier  was 
determined  by  taking  the  ratio  of  the  scattered  radiation  dose  values  with  the 
horizontal  barrier  in  position  to  the  respective  measured  values  for  the  zero  psf 
horizontal  barrier  case.  Measured  attenuation  values  versus  average  solid  angle 
fraction,  w,  are  shown  in  Figure  5.  This  average  solid  angle  fraction,  u,for  the 
vertical  wall  subjected  to  the  20.5  in.  radius  incident  beam  was  approximately 
determined  from  u|  and  <Jg  values  based  on  the  lower  and  upper  limits  of  an  equivalent 
square  area  (36  by  36  in.)  of  incident  radiation.  Figure  5  also  presents  the  theoretical 
attenuation  curves  for  36.5  and  75  psf  ceilings,  as  calculated  by  the  method  used  in 
this  report  for  cesium  (.66  mev)  radiation  impinging  perpendicularly  on  the  vertical 
slab.  The  majority  of  the  experimental  points  fall  below  the  theoretical  curves 
except  at  small  values  (w-  .175)  of  <*J,  where  experimental  dose  rates  were  higher 
than  those  predicted  by  theory.  This  increase  above  the  theoreticol  curves  at  small 
u  is  not  expected  and  cannot  be  explained  at  this  time.  The  values  of  dose  upon 
which  these  data  at  small  u  are  based  are  extremely  low  and  any  unaccounted  air- 
scattered  or  ground-penetration  components  even  in  small  quantities  could  contribute 
significantly  at  these  low  rate  regions,  providing  a  false  reading  of  this  type.  Results 
of  other  experiments  discussed  later  in  this  chapter  did  not  show  this  crossover  at 
small  values  of  <*. 


"In  and  Down"  Measurements  at  the  Protection  Structures  Development  Center 

Cobalt-60  measurements  were  performed  on  a  24  by  36  ft.  rectangular 
concrete  structure  with  a  basement  having  an  exterior-wal I  thickness  of  50  psf. 

The  experiment  was  conducted  both  with  and  without  a  50  psf  concrete  basement 
ceiling  slab,  so  that  the  basement  ceiling  effect  could  be  measured  directly.  In 
Figure  6  is  presented  the  comparisor.  between  the  measured  ceiling  effect,  Bo(Xf,u) 
and  that  taken  from  the  theoretical  values  of  this  report  as  given  in  Figure  3.  The 
50  psf  theoretical  curve  falls  slightly  above  the  50  psf  measured  values,  giving  a 
tatio  of  experiment  to  theory  of  approximately  0.76.  The  computed  values  in  this  case 
therefore  were  both  conservative  and  close  to  the  0.815  mean  experiment  to  theory 
ratio  determined  for  c'l  the  experiments,  covered  in  this  report.  (See  concluding 
section  of  this  chapter.) 


Attenuation  of  Cobalt-60  Radiation  by  a  Small  Iron  Cylindrical  Structure  ' 

A  series  of  measurements  were  conducted  on  a  2  ft.  diameter  cylindrical 
iron  structure.  The  structure  was  composed  essentially  of  an  iron  pipe  2  ft.  in 
diameter  mounted  vertically  in  the  ground  so  that  only  the  upper  2  ft.  of  its  length 
projected  above  the  ground  level.  The  above-ground  wall  thickness  of  this  structure 
was  varied  from  1/4  to  1-1/2  in.  (approximately  5-60  psf)  with  floor  thicknesses  of 
0,  10,  20,  and  40  psf  iron.  Below-ground,  i.e.,  basement  measurements  were  made 
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along  the  vertical  centerline  at  depths  of  1/4,  1/2,  3/4,  1,  2,  and  3  ft.  Results 
expressed  as  the  ratio  of  the  measured  dose  rate  in  the  basement  with  the  ceiling  in 
place  to  that  measured  without  the  ceiling  (zero  psf),  a  direct  measure  of  attenuation 
factor  Bo(Xf,u),  are  presented  in  Figure  7. 

The  measured  floor  reduction  factors  vary  with  solid  ongle  fraction  in  a 
manner  similar  to  that  of  the  theoretical  curves  based  on  the  theory  developed  in 
this  report.  The  measured  values  are  somewhat  lower  (more  conservative)  than 
predicted  by  theory,  as  would  be  expected,  since  all  assumptions  and  approximations 
used  in  developing  the  theory  were  aimed  at  producing  a  conservative  attenuation 
value. 

Basement  Attenuation  —  Other  Experiments 

A  comparison  of  experimentally  measured  dose  rates  with  those  computed 
by  the  Engineering  Manual^  method  (modified  for  Cobalt-60  concrete),  using  both 
the  floor  attenuation  data  given  in  that  report  and  those  provided  by  the  method 
developed  in  this  report,  is  presented  in  Table  4  for  all  currently  reported  experiments. 

In  general,  agreement  between  expei (mentally  measured  values  and  the  treatment 
recommended  by  this  study  is  much  better  than  that  provided  by  the  method  currently  used. 
In  particular,  the  basement  ceiling  barrier  experiments  currently  being  made  at  NDlJO 
(on  a  10  by  20  ft.  roofless  concrete  structure  with  46  psf  exterior  walls  and  floor  slab), 
though  still  preliminary,  provide  much  better  agreement  with  the  method  used  in 
this  report  than  Engineering  Manual  values,  particularly  at  the  lower  depths.  The 
same  trend  holds  true  for  the  data  in  an  older  study^  ^  of  a  simple  structure  with  a 
basement,  even  though  those  data  are  somewhat  difficult  to  interpret  because  of 
inhomogeneity  in  floor  mass.  The  mean  ratio  of  the  experimental  results  of  these 
two  full-scale  studies  to  those  using  the  computed  method  of  this  report  is  approximately 
0.70,  compared  to  a  mean  of  0.815  for  all  experiments  reported.  Values  calculated 
by  the  method  used  in  this  report  are  conservative  compared  with  those  provided  by 
experiment  and  follow  the  dose  versus  depth  contours  better  than  previous  analyses. 

Protection  factor  versus  depth  contour  agreement  is  also  good  for  the 
Kansas  State  University^  ^  blockhouse  end  the  Technical  Operations  multistory  1/12 
scale  steel  modelsJ  '  ^  The  KSU  structure  was  a  square  20  by  20  ft.  concrete 
house  of  60  psf  walls  and  56  psf  floor  and  roof;  the  T/O  model  was  of  a  36  by  48  ft. 
six-story  windowless  building.  The  absolute  measured  level  of  the  dose  in  the 
basement  of  the  KSU  blockhouse  is  higher  (prediction  less  protection)  than  that 
computed  by  the  method  of  this  report  -the  only  time  that  this  has  occurred  in  a 
full-scale  structure.  Dose  rates  in  the  basement  of  the  80  psf  wall  and  floor  steel  model 
(T/O)  also  were  higher  than  predicted.  Reported  misalignment  of  the  structure  with  the 
basement  and  edge  scattering  of  radiation  from  the  relatively  thick  basement  ceiling 
s'^b  that  projected  above  the  ground  level  indicates,  however,  that  the  absolute  value 
of  !  e  d~> se  measured  is  of  questionable  validity  in  this  case. 
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!  TABLE  4 

(Dose  Rate  in  Belov/  Ground  Regions  from  Ground  Based 

Sources  of  Contamination 


>epth,  ft 

Experiment 

Engineering 

Manual 

Method  of 

This  Report 

NDL10 

1 

.0034 

.0035 

.0052 

2 

.0033 

.0023 

3 

.0016 

4 

.0029 

.0010 

.0043 

5 

.0026 

.00068 

.0035 

6 

.0024 

. 00048 

.0028 

Simple  structure 

1-3/4 

.0019 

.0010 

.0025 

with  basement^  ' 

4-1/4 

.0018 

. 00098 

.0023 

6-3/4 

.0017 

.00101 

.0026 

9-1/4 

.0015 

.00040 

.0020 

12 

KSU  blockhouse 

1/4 

.0041 

.0026 

.0025 

1-3/4 

.0048 

.0018 

.0025 

Multistory  steel 
m„deli3'H 

4-3/4 

.0051 

. 00082 

.0026 

20  psf  walls 

3 

.0195 

.0142 

.0259 

20  psf  floors 

10 

.0185 

.0055 

.0185 

15 

.0140 

.0035 

.0106 

20  psf  walls 

1 

.00094 

.00164 

.00101 

80  psf  floors 

3 

.00115 

.00133 

7 

.0016 

.00075 

10 

.0018 

.00052 

.00176 

15 

.0018 

. 00034 

.000127 

80  psf  walls 

1 

.00085 

.00072 

.00041 

80  psf  floors 

3 

.0030 

. 00058 

.00051 

7 

.0033 

. 00034 

. 00055 

10 

.0034 

. 00023 

. 00070 

15 

.0026 

.00015 

.00056 

40  psf  wails 

3 

.0022 

.0034 

.0040 

50  psf  floor: 

6 

.0023 

i  . 0022 

.0040 

9 

.0021 

.0017 

.0048 

15 

.0022 

. 00087 

.0033 

t 
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Figure  7 
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Summary  of  Comparisons  with  Experiment 


The  method  o:  this  report  for  computing  basement  protection  factors,  using 
the  value  Bo(Xf,w)  in  determining  the  effect  of  the  basement  ceiling  barrier,  provides 
results  that  agree  much  better  with  measuit’d  basement  attenuation  than  does  the 
Engineering  Manual  method,  particularly  as  detector  depth  below  the  ceiling  slab 
increases.  When  comparing  computed  values  based  on  the  method  used  here  with 
measured  values,  nearly  all  the  computed  values  indicated  slightly  less  attenuation 
than  the  comparable  measured  values.  The  computed  or  theoretical  values  therefore 
are  conservative  as  compared  with  measured  values. 

Figure  8  illustrates  the  comparison  between  measured  and  computed  values 
of  Bo(Xf,u)  of  this  reporr  for  the  three  experiments  providing  data  for  basement 
positions  both  with  and  without  the  ceiling  slab  in  place  (i.e.,  a  direct  measurement 
of  barrier  factor  B'0).  All  values  fall  on  the  conservative  side  of  the  figure  except 
for  a  few  of  the  DCBRL  points  token  at  small  solid  angle  fractions  (u  <.175),  where 
the  effects  of  any  radiation  component  not  accounted  for  (skyshine  is  suspected)  in 
these  low  dose  rate  areas  would  provide  a  significant  part  of  the  measured  dose  and 
would  therefore  seem  to  reflect  a  major  decrease  in  the  floor  slab's  effectiveness. 
Figure  9  shows  the  distribution  of  the  ratios  of  all  measured  basement  attenuation 
experimental  values  compared  with  that  computed  using  the  theory  of  this  report. 

At  ratios  greater  than  1.10,  the  tail  on  this  graph  can  be  attributed  almostentirely 
to  measurements  at  small  solid  angle  fiactions  for  the  DCBRL  experiment  mentioned 
above,  to  the  KSU  data,  and  to  the  questionable  80  psf  wall  and  floor  T/O  model 
multistory  structure  data.  A  statistical  analysis  of  this  distribution  is  shown  in 
Figure  10,  which  clearly  indicates  the  mean  value  of  the  experiment  to  computed 
dose  ratio  is  0.815  with  a  sfcndard  deviation,  a,  of  0.135  for  all  the  experimental 
data. 


Since  Figure  10  represents  a  sizable  number  of  experimental  data 
(approximately  350  experimental  points),  it  is  recommended  that  the  mean  ratio  of 
experiment  to  theory  of  0.815  be  used  to  reduct  the  computed  barrier  function 
(Figures,  2,  3,  and  4)  of  this  report  in  order  to  provide  better  agreement.  Clearly, 
it  is  not  enough  to  multiply  tire  computed  values  by  this  ratio,  because  the  attenuation 
of  a  zero  thickness  floor  must  remain  1.0.  In  Figure  17  of  Chapter  3,  the  floor 
barrier  attenuation  for  the  1.12  hr.  fallout  spectrum  has  been  corrected  by  multiplying 
by  0.815  for  floor  thickness  greater  than  20  psf,  by  .91  for  floor  thickness  of  10  psf 
and  faired  into  a  vo'ue  or  1.0  for  floor  thicknesses  of  zero  psf.  This  chart  thereby  is 
brought  into  better  agreement  with  expected  experimental  values.  It  is  further 
recommended  that  a  similar  procedure  be  used  for  cobalt  and  cesium. 

A  report  entitled  "Dose  Albedo  and  Transmission  Coefficients  for  Cobalt-60 
and  Cesium-137  Gamma  Rays  Incident  on  Concrete  Slabs'1^  by  M.  J.  Burger  and 
E.  E.  Morris  was  published  when  this  report  was  nearing  completion.  The  attenuatior, 


21 


Theoretical  Attenuation  Factor  B‘  (Xr,  u  ) 


Measured  Attenuation  Factor 


Theory  Versus  Experiment 


Figure  8 
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of  parallel  slant  incident  radiation  for  the  values  used  in  calculating  the  floor 
attenuation  for  this  work  at  1.25  mev  is  compared  with  *hot  of  Spencer^  ct  1.25  rr.ev 
and  Raso  at  1.00  and  2.00  mev  in  Figure  11  for  o  typical  case.  This  figure 
presents  the  comparable  attenuation  curves  as  a  function  of  the  cosine  of  the  slant 
incidence  angle  0  for  a  slab  of  one  mean  free  path  thickness.  This  latest  work  of 
Berger  and  AAorris  yields  lower  values  of  attenuation  (i.e.,  the  transmitted 
radiation  level  is  lower).  The  Berger  end  Moms  attenuation  coefficients  ,in  general, 
are  lower  by  10  to  20  per  cent  than  an  interpolated  1.25  mev  cur/e  for  Raso’s  work, 
the  greatest  difference  occurring  in  the  range  0.3  <  cos  0  <0.6.  The  attenuation 
factor  Bg(Xf,  u)  proposed  in  this  report  for  use  in  basement  calculations  was  based 
on  the  transmission  data  of  Raso.  It  was  further  proposed  that  the  Bo(Xf,«)  results 
be  reduced  by  o  factor  of  0-815  to  bring  theory  and  experiment  into  closer  agreement. 
If  Berger  and  Morris  transmission  data  hod  been  used  in  place  of  Raso's,  the 
theoretical  value  of  Bo(Xf,u)  would  have  agreed  nuch  more  closely  with  experiment 
without  the  0.615  factor.  This  confirms  the  necessity  of  the  0.815  reduction  factor 
to  bring  the  analysis  in  line  with  experiment. 


25 


C  O  N  ESC  O  consultants  in  engineering  science 


Supplement  to  TR-20-Vol.  2 


CHAPTER  3 

DRAFT  SUPPLEMENT  TO  TR-20  Vol.  2  and  PM100-1 
FOR  IN  AND  DOWN  CORRECTION 

IMPROVED  BASEMENT  CALCULATION 


INTRODUCTION 


Recent  studies  have  indicated  that  more  accurate  protection  factors  may 
be  achieved  in  basement  areas  if  the  attenuation  of  the  basement  ceiling  is  con¬ 
sidered  to  be  a  function  of  both  ceiling  thickness  and  structure  size.  In  this  memorandum 
we  shall  summarize  briefly  the  method  of  calculating  the  attenuation  provided  by  a 
basement  ceiling  and  will  provide  the  charts  and  graphs  necessary  to  make  this 
calculation. 


Floor  Barrier  Corrections,  &XJ(Xf  ) 

In  the  current  edition  of  the  equivalent  building  method,  the  floor  barrier 
effect  is  represented  only  as  a  function  of  the  floor  thickness.  Recent  experiments 
indicate  that,  in  some  cases,  this  interpretation  may  be  somewhat  in  error.  New  barrier 
functions  providing  much  better  agreement  have  been  computed  by  assuming  single 
Compton  scatter  energy  loss  in  the  structure  walls  and  attenuation  similar  to  that 
experienced  by  parallel  slant  incident  photons  through  the  floor  slab.  Floor  barrier 
effect  thus  becomes  a  function  of  both  floor  thickness  and  structure  area. 

Figure  ll-l  of  TR-20  Vol.  2  indicates  that  the  basement  Pf  charts  are  based 
on  a  floor  barrier  of  1.0  or,  equivalently,  a  mass  thickness  of  zero  for  the  floor  above 
the  detector  This  was  done  so  that  any  corrections  would  be  additive.  The  effect 
of  basement  ceiling  thickness  in  attenuating  radiation  fiom  ground-based  sources  of 
contamination  is  presented  in  auxiliary  charts,  one  for  each  basement  area. 

New  figures  replacing  Charts  5  through  9  of  TR-20,  Vol.  2  accompany  this 
supplement.  When  these  charts  are  used  to  compute  Pf  for  a  bsement  location,  the 
auxilary  charts  on  the  upper  poition  of  the  figure  must  be  used  to  correct  for  the 
barrier  effect  of  the  floor  above  the  detector.  The  auxiliary  charts  were  derived  by 
converting  for  each  building  area  the  attenuation  values  computed  as  described  above 
to  equivalent  weight  of  wall  barrier  factor  Be  from  Chart  2  of  the  Engineering  Manual 
(TR-20,  Vol.  1). 
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Illustrative  Example 

What  is  the  protection  factor  in  the 
basement  for  a  rectangular  building 
90  ft.  wide  by  1 10  ft.  length. 

Floor  thickness  Xf  =  50  psf 

Wall  thickness  Xw  =  70  psf 

Roof  thi  ckness  Xr  =  100  psf 

Solution: 


X 


Building  area  =  90  x  110'  =  9900  2  ~  10,000  sq.  ft.;  therefore,  Chart  8 

/1 0 1 

Equivalent  basement  area  Ag  =  9900  =  3,430  sq.  ft. 

1  o2 

Equivalent  roof  area  A'  =  9900  =  2250  sq.  ft. 


Therefore,  from  Chart  16  (TR-20  -  Vol.  2),  equivalent  roof  thickness  ~  152  psf 
Wall  equivalent  ceiling  thickness  AXU  (Xf)  =  195;  Chart  8,  top 
Enter  Chart  8  with  xw  =  70  +  195  =  265  psf  and  XQ  =  152 

Read  Protection  Factor,  Pf  =  180 
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Supplement  to  PM  100-1 
TR  20  Vol.  1 


DRAFT 

IMPROVED  "IN  AND  DOWN"  CORRECTION 


INTRODUCTION 


Recent  studies  have  indicated  that  protection  factors  of  higher  accuracy 
may  be  achieved  in  basement  areas  if  the  attenuation  of  the  basement  ceiling  is 
considered  to  be  a  function  of  both  ceiling  thickness  and  structure  size.  In  this 
memorandum  v/e  shall  summarize  briefly  the  method  of  calculating  the  attenuation 
provided  by  a  ceiling  and  will  provide  a  new  chart  of  floor  barrier  factor  for  ground 
contributions  through  the  floor  above  the  detector  to  make  this  calculation  possible. 


Floor  BarrierFactor,  Bq 

In  the  current  edition  of  the  professional  manual  entitled  "Shelter  Design 
and  Analysis,  TR-20  -  Volume  1,  Fallout  Protection  and  PM-100-1,  Design  and 
Review  of  Structures  for  Protection  from  Fallout  Gamma  Radiation,"  the  attenuating 
effect  of  the  floor  slab  is  represented  as  a  function  of  only  the  floor  thickness. 

Recent  experiments  indicate  that  in  some  cases  this  interpretation  may  be  somewhat 
in  error.  New  barrier  functions  providing  much  better  agreement  have  been  computed 
by  assuming  single  Compton  scatter  energy  loss  in  the  structure  walls  and  attenuation 
similar  to  that  experienced  by  parallel  slant  incidence  photons  through  the  floor 
slab.  Floor  attenuation  Bq  thus  becomes  a  function  of  both  floor  thickness  and 
structure  geometry. 

This  new  attenuating  function,  Bo(Xq,w),  replacing  Bq(X'0),  is  used  in 
the  same  way  as  that  previously  used.  The  expected  dose  rate  in  a  structure  from 
wal  I -scattered  sources  of  radiation  is  computed  for  a  floor  slab  of  zero  thickness 
and  then  multiplied  by  the  attenuation  introduced  by  the  floor.  Since  the  attenuation 
afforded  by  the  floor  is  a  function  of  the  solid  angle  fraction  of  the  scattering  mass 
relative  to  the  protected  position  of  interest,  as  well  as  the  floor  thickness  (see 
the  figure  accompanying  the  illustrative  example),  an  exact  calculation  would 
require  averaging  over  the  total  scattering  wal!  from  Uy  to  An  approximate 
method  that  provides  excellent  agreement  with  this  more  exact  procedure  is  to 
determine  and  use  a  value  of  attenuation  Bo(Xo/u)  for  the  arithmetic  average  of 
the  solid  angle  fraction  describing  the  scattering  surface. 
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As  in  the  illustrative  example,  the  attenuation  introduced  by  the  floor  would 
be  equivalent  to  B^(X'q,  u  ),  where: 


1.  Functional  equations 

Cq  =  negligible  because  XQ  is  greater  than  350  psf 
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2.  Determining  solid  angle  fractions  and  parameters 


Chart 

3* 

Chart 

5* 

Chart 

5* 

W 

r~ 

z 

e 

u 

GaM 

G(u) 

u 

u 

90 

no 

24 

.82 

. 

.44 

.33 

w' 

u 

90 

110 

10 

.82 

.18 

.82 

m 

u  +  u* 

—  _  u  u  .60  +  .82 

Y~  ^ - 


B'  (X  ,u 
o  c 

=  8^(50,0.71) 

=  0.029 

Accompanying  chart  (Fig.  17) 

B  (X  h) 
e  e 

=  Be(70,3) 

=  0.18 

Chart  No.  2* 

S  (X  ) 

w  e 

=  S  (70) 

w 

=  0.66 

Chart  No.  7* 

E(e) 

-  £(.82) 

=  1.41 

Chart  No.  8* 

♦Chart  of  PM  100- 1 


3.  Solution: 


C  =  (0.18)  (0.02?)  {  T (.079  -  .051)  (1 -.66)  + 


.33-  .19 


.661 


1.41 


C  =  .00073 

g 

K  =  c  +  C  =  .00073  + .00000  =  .00073 
r  g  o 

Pp  =  l/Rp  =  1/ .000 73  =  1370  say  1300  answer 
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Illustrative  Example  2  -  Ref.  Example  5.3, 
TR-20  (Vol.  1) 

Exposed  Basement  Walls 


L 

100  ft. 

W  =  50  ft. 

X 

80  psf 

(first  story) 

e 

X  '  = 

100  psf 

(basement) 

e 

x»  = 

60  psf 

0 

X 

100  psf 

r 
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Chart 

3 
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5 

Chart 

5 

u 

w 

L 

z 

e 

n 

U) 

G 

s 

G 

a 

u” 

u 

50 

100 

23 

.50 

.46 

.47 

.39 

.036 

u1 

u 

50 
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11 
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Soluti  on: 


a. 


u 


u‘  +u" 

u  u 
- 


.70  +  .47 

2 


.59 


B  (X  ,H)  =  .095 
e  e 

S  (X  )  =  .73 

w  e 

E(e)  =1.34 


C 

9 


G  (u1  ) 
a  u 


G  (u  ) 
a  g 


S  (X  ) 

w  e 


+ 


G  (o'  )  -  G  (u  ) 
s  u  s  u 


i  (X  )E(e)>  B  (X  ,H) 

we  i  e  e 


C 

9 


(.069  -  .029)  (1-.73)  +  (.28  -.11)  (.73)0.34) 


C  =.  .0168 
9 


.095 
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c.  Roof  Contribution 

v  =  X  +  X’  =  160  psf 

o  r  o 

i 

C  (u\X  )  =  .0039 

o  U  0 

Pf  =  47 

Detector  heights  above  ground  are  treated  in  a  similar  manner 
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CHAPTER  4 

FINITE  FIELDS  OF  CONTAMINATION 


INTRODUCTION 


Many  shelter  areas  in  this  country  depend  in  part  for  their  worth  on  the 
fact  that  they  are  subjected  to  limited  rather  than  infinite  fields  of  contamination. 

A  common  situation,  that  of  row  housing  in  urban  areas,  for  example,  would 
subject  the  sheltering  structures  to  radiation  originating  from  a  finite  field  of 
contamination  perhaps  50  ft.  wide  on  both  sides  of  the  structure.  These  structures 
would  then  be  as  much  as  three  times  more  effective  than  an  identical  structure 
surrounded  by  a  full  field  of  contamination.  Parker^  estimates  that  fully  half 
of  the  basement  shelter  areas  surveyed  at  the  time  his  report  was  published  were 
subjected  to  field  widths  of  less  than  60  ft. 

Two  areas  of  inquiry  arise  when  we  evaluate  the  methods  currently  used 
in  predicting  the  effects  of  finite  fields  of  contamination.  We  must  first  determine 
if  the  more  detailed  method  of  calculation  as  represented  by  the  "engineering 
method"  adequately  portrays  the  physical  situation.  Second,  we  must  ask  whether 
or  not  the  simplified  "equivalent  building  method"  agrees  with  results  obtained 
by  means  of  the  more  complex  method.  The  agreement  between  the  two  methods  in 
their  present  form  has  been  examined^  previously  and  recommendations  for  improving 
upon  this  agreement  have  been  made.  The  required  changes  have  been  held  in 
abeyance  until  an  adequate  answer  to  the  first  area  of  inquiry  (the  accuracy  of 
the  more  complex  method  of  calculation)  has  been  determined.  Our  purpose  in 
this  chapter  is  to  examine  this  accuracy  by  comparing  if  with  the  best  information 
available. 

There  are  few  experimental  data  with  which  to  evaluate  the  methods 
proposed  for  computing  the  shelter  afforded  against  finite  fields.  Several  years 
ago,  an  experiment  was  conducted  to  determine  the  effect  of  finite  fields  of 
contamination  on  above-ground  locations.  This  experiment  used  the  modeling 
technique,  with  iron  as  the  attenuating  medium,  and  simulated  contamination  in  the 
form  of  rectangular  areas.  Eisenhauer*  analyzed  these  data  in  detail,  concluding 
that  "comparison  of  experimental  and  calculated  results  shows  good  agreement  for 
both  structures  except  for  the  first  story."  A  second,  brief  experimental  series,  con¬ 
ducted  on  full-scale  concrete  structures  by  McDonnell,  using  circular  annuli, 
again  indicated  fair  agreement  between  experiment  and  theory.  Both  experiments 
were  exclusively  for  upper-story  locations,  where  both  direct  and  scattered  radiation 
affected  the  dose  measured  in  approximately  equal  amounts.  Since  neither  of  these 
experimental  series  was  rigorous  in  the  examination  of  finite  fields  of  contamination 
and  because  so  many  shelter  areas  seem  to  be  affected  only  by  finite  fields  of 
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contamination,  a  further  evaluation  was  undertaken.  Eisenhauer  recently  com¬ 
pleted  a  computer  compilation  program,  using  Berger  and  Morris'  Monte  Carlo 
transmission  coefficients '  ^  for  1.25  mev  prediction.  This  program  is  capable  of 
computing  the  raaiction  originating  from  any  arbitrarily  located  square  patch  of 
contamination  through  any  rectangular  wall  area  to  an  arbitrarily  located  detector. 
The  program,  therefore  may  be  used  to  generate  numerical  "experiments"  with 
which  to  compare  the  existing  calculational  method. 


The  Engineering  Method  of  Shelter  Computations  from  Finite  Contamination  Fields 


The  engineering  method  of  shelter  computation  described  in  the  manual 
"Design  and  Review  of  Structures  for  Protection  from  Fallout  Gamma  Radiation"^ 
presents  rules  for  computing  the  ground-based  dose  contributions  for  structures 
subjected  to  radiation  from  both  infinite  and  limited  fields  of  contamination.  In  this 
method  the  radiation  arriving  at  a  point  within  a  structure  is  subdivided  into  three 
components:  radiation  that  (1)  has  passed  directly  through  the  building  walls  without 
scattering,  (2)  has  been  scattered  by  the  wall',  and  (3)  has  been  scattered  by  the 
atmosphere.  Non-wall  scattered  or  direct  radiation  from  an  infinite  field  of  con¬ 
tamination  is  determined  by  multiplying  the  cumulative  angular  distribution  of 
non-wall  scattered  radiation,  Gj^h)1  as  viewed  from  the  point  of  interest  in 
the  structure  by  the  height-dependent  wall  barrier  factor  B(Xe,h)i  and  by  the 
fraction  of  ladiation  not  scattered  by  the  structure  walls,  (1-SW).  Calculations 
for  the  finite  field  case  are  similar  except  that  the  cumulative  angular  distribution 
of  non-wall  scattered  radiation  must  be  differenced  to  account  for  the  finite  field. 

The  wal  I  -scattered  component  for  an  infinite  field  case  is  determined  by 
multiplying  the  cumulative  angular  distribution  of  radiation  scattered  from  the 
structure  walls,  Gs(u)  ,  by  the  height-dependent  wall  barrier  factors  B(Xe,h) 
the  shape  correction  facror,  (E),  the  fraction  of  emergent  radiation  scattered,  (Sw), 
and  a  floor  barrier  factor,  B0(Xf )  or  ceiling  factor,  Bo(Xc  )  ,  where  applicable. 

The  floor  barrier  factor  is  used  for  attenuating  wal  I -scattered  radiation  reaching 
the  detector  from  the  floor  below  in  a  multistory  structure  whereas  similarly,  the 
ceiling  barrier  factor  is  used  for  wall-scattered  radiation  reaching  the  detector 
from  the  story  above  For  a  finite  field  case  this  procedure  is  modified  by  applying 
a  different  wall  barrier  factor,  Bu^u^Xg)  ,  which  is  expressed  as  a  function  of 
the  solid  angle  of  the  field  of  contamiation  as  viewed  from  the  structure's  wall  at 
midsto'y  height  in  place  of  the  infinite-field  wall-barrier  factor,  B(Xe,h). 

The  atmospheric  scattered  component  for  both  the  finite  and  the  infinite 
field  of  contamination  cases  are  determined  by  multiplying  the  cumulative  angular 
distribution  of  skyshine  radiation,  Ga(w),  ,  by  the  wall  barrier,  B(Xe,h)  ,  the 
fraction  of  radiation  not  scattered  in  passing  through  the  structure  wails,  (1-SW), 
and  the  ceiling  barrier,  Bo(Xe)  ,  when  appropriate.  The  assumption  that  the 
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atmospheric  scattered  component  for  a  finite  field  is  nearly  identical  to  that  of 
the  infinite  field  case  should,  in  practice,  be  fairly  accurate.  Though  adjacent 
buildings  may  block  the  direct  view  of  part  of  the  field  of  contamination,  however, 
the  air-scattered  radiation  physically  blocked  by  the  adjacent  buildings  will  be 
approximately  balanced  by  radiation  scattering  from  the  walls  of  these  buildings. 
Also,  it  is  expected  that  in  most  realistic  cases  .here  will  be  fallout  on  the  tops 
of  the  adjacent  buildings  as  well  as  in  the  region  behind  them;  hence  the 
angular  distribution  of  the  radiation  that  is  not  blocked  by  adjacent  walls  will 
remain  the  same. 


The  three  radiation  components:  direct,  wail -scattered,  and  air-scattered, 
are  then  summed  to  give  the  total  dose  contribution.  The  general  equation  for  the 
total  ground  contamination  within  a  multistory  structure  in  a  limited  field  of 
contamination  is  presented,  together  with  a  sketch  of  the  idealized  building 
arrangement,  in  Figure  18. 


Charts  in  the  Engineering  Manual  for  determining  the  geometric  and 
barrier  reduction  factors  used  in  structure  shielding  calculations  were  derived  from 
the  basic  data  on  radiation  penetration  developed  by  Spencer'^  for  1 .  ] 2  hr.  fallout. 
Chart  7  in  the  manual,  for  the  fraction  of  emergent  radiation  from  a  wall  barrier 
that  is  scattered,  (Sw),  is,  however,  for  Cobalt-60  radiation.  Since  most  experi¬ 
mental  verifications  of  the  Engineering  Manual  calculation  procedures  have  been 
and  are  being  performed  using  Cobalt-60  radiation,  and  since  the  NBS  computer 
compilation  program^  against  which  we  wish  to  compare  the  engineering  method 
uses  the  1.25  mev  (Cobalt-60)  Monte  Carlo  transmission  coefficients  of  Berger 
and  Morris,  the  Engineering  Manual  method  has  had  to  be  extended  to  1.25  mev 
radiation  in  order  to  permit  accurate  comparisons.  Cobalt-60  curves  for  the 
barrier  shielding  reduction  factors,  Bo(Xf)  Bo(Xe,3‘),  and  Bq(Xc)  are  presented 
in  Figure  19a.  The  data  in  this  figure  has  been  calculated  in  the  same  way  and 
are  equivalent  to  Chart  1  of  the  Engineering  Manual,  in  which  Case  1  applies  to 
a  horizontal  floor  barrier,  Bo(Xf),  Case  2,  the  exterior  vertical  wail  barrier, 

B(Xe,3'),  and  Case  3,  an  overhead  or  ceiling  horizontal  barrier,  Bq(Xc).  Figure  19b 
presents  the  Cobalt-60  overhead  or  ceiling  horizontal  barrier  Bo(Xc).  Figure  19b 
presents  the  Cobalt-60  exterior  wall  attenuation  effects  for  detector  heights  to 
300  ft.,  and  is  equivalent  to  the  Chart  2  fallout  curves  of  the  Engineering  Manual. 
Cobalt-60  directional  responses  from  ground  sources  of  contamination  Gj,  Gs, 
and  Ga  are  given  in  Figure  20a,  which  has  been  calculated  in  the  same  way  and 
therefore  is  equivalent  to  Chart  5  of  the  Manual  for  1.25  mev  radiation.  Figure  20b 
similarly  extends  the  directional  response  for  direct  radiation  G(J  to  detector  heights 
of  up  to  100  ft.,  and  is  equivalent  to  Chart  6  of  the  Manual. 


Barrier  shielding  effects  B^Xf),  Be(Xe,  and  Bq(Xc)  for  Cobalt-60,  os 
given  in  Figures  19a  and  19b,  were  developed  directly  from  Spencer's^*  L(X), 
2W(X,d),  and  S'(X)  curves  respectively,  and  agree  closely  with  those  developed  by 
LeDoux.^l 
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Figures  20o  and  20b  for  directional  response  Gj(u,h),  Gs(u)  ond  GQ(X) 
were  developed  from  the  basic  work  of  Spencer.  *6  LeDoux^l  also  has  developed 
Cobalt-60  curves  for  the  directional  responses,  which  in  general  agree  well  with 
the  curves  of  Figure  20  except  for  values  of  G<j  at  detector  heights  above  3  ft. 

Values  for  the  direct  angular  distribution  G^(u,h)  in  Figures  20a  and 
20b  were  developed  from  the  equation: 


cos  9v 

r 

j  (d,cos9)d(cos0) 

G>,h)=  9-, - 

I  K  (d,co s9)d(cos9) 


1  o  .1 

(d,cos9)d(cos9)  -  J<(d,co$0)d(cos9)  -J^(d,cos0)d(cos0) 
- 1  - 1  co$9x, 

^,,,(d/cos9)d(cos9) 


which  in  Spencer's^  notation  is  identically  equal  to 


Gd(u,h)  =  1 


5(d) 

Lid) 


-  L  (d,u) 
a 


where  t'  (d,  cos  9)  is  the  angular  distribution  and  S(d),  L(d),  and  La(d,u)  for  Cobalt-60 
are  taken  from  Spencer's^  work.  S(d)  divided  by  L(d)  is  the  ratio  of  skyshine  to 
total  (adiation  at  a  distance  d  above  a  uniform  infinite  source  field.  Gs(cobalt) 
was  calculated  from  the  equation  used  for  Gs(fallout),  using  Sa(d,u)  curve  for 
Cobalt-60  at  d  =  3  ft.  Thus: 


G  =  0.5  (1-S  (d,u) 
s  a 


Similarly,  the  cumulative  skyshine  contribution 


G  =  0.168  G  (u)  1  +  (1/2  -  G  (u) 


where  0.168  is  twice  the  skyshine  fraction  of  the  total  radiation  at  3  ft.  altitude. 
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Agreement  for  Gjfu^1)  between  values  calculated  in  this  manner  and 
curves  developed  by  LeDoux  are  excellent  up  to  u  ;  .90.  For  values  of  o  -.90, 
the  values  of  Figure  20  are  somewhat  higher  than  those  of  LeDoux  (for  example, 
at  u  -  .99,  Gd  (Figure  20)  .14,  G^(LeDoux)  -  .105).  Similarly,  agreement  for 

the  cumulative  angular  distribution  of  wall -scattered  radiation  G$  between 
Figure  20  and  values  previously  derived  by  LeDoux  is  excellent,  Conesco  values 
being  slightly  higher  at  o  .75  (for  example,  at  u  =  .99,  Gs(Figure  20)  =  .013,  Gs 
(LeDoux)  =  .01  1).  The  skyshine  cumulo  ve  angular  distribution  Ga  also  shows 
good  agreement  between  the  values  shown  in  Figure  20  and  those  previously 
derived  by  LeDoux,  with  values  of  Figure  20  slightly  less  for  u  .75  and  slightly 
higher  for  u  .75.  Figure  20b  presents  the  direct  response  Gj  for  heights  to 
100  ft.  and  agrees  fairly  well  with  values  developed  by  LeCoux.  Agreement  at 
Gd(u,  3')  as  discussed  above,  was  good.  Agreement  at  heights  greater  than  3  ft, 
however,  is  not  as  good,  particularly  at  small  values  of  o.  This  is  illustrated  in 
Figure  20  by  including  both  curves  at  Gd(u,h)  equal  to  0.50  and  0.75. 

Engineering  Manual  Chart  9,  used  for  obtaining  the  finite  field  wall 
barrier  factor  Bus  (us,Xe)  is  for  1.12  hr,  fallout.  Fallout  values  for  Bus  (Us/Xe) 
were  used  in  this  study  for  calculating  1.25  mev  structure  reduction  factors  for 
comparison  with  the  computer  Monte  Carlo  results.  This  was  done  for  two  reasons: 

(1)  a  Cobalt-60  barrier  chart  for  limited  strips  of  contamination  had  not  been 
developed  at  the  time  this  study  was  conducted;  (2)  the  error  introduced  by  using 
the  fallout  chart  values  in  place  of  Cobalt-60  should  be  small  at  the  mass  thickness 
(9,  36,  and  144  psf)  used  in  this  study.  The  error  in  the  finite  wall-barrier  factor 
in  using  failoui  values  instead  of  Cobalt-60  can  be  approximated  by  considering 
the  differences  between  fallout  and  Cobait-60  infinite  field  wall  barier  factors. 

For  the  infinite  field  case  the  Cobalr-60  barrier  factor  is  about  5  per  cent  less  than 
the  fallout  value  at  100  psf,  and  25  per  cent  at  200  psf.  Since  the  wall  thickness 
values  used  in  tfvs  study  are  36  psf  or  less,  except  for  one  point  at  144  psf,  errors 
introduced  are  expected  to  be  small. 

Computations  of  the  Engineering  Manual  type  were  performed  using  data 
for  1.25  mev  on  a  series  of  simple  structures  for  comparison  with  corresponding 
values  determined  by  the  NBS  compilation  prog_om  for  finite  fields.  T.ie  series 
v/as  based  on  a  single-story  structure  of  the  blockhouse  type.  Building  floor  area, 
widtn-to-length  ratios  (eccentricity),  and  wall  thickness  were  varied  to  provide 
data  for  a  comparison  with  the  compilation  program  over  a  wide  rangr*  of  structure 
parameters.  The  width  of  the  rectangular  field  of  contamination  surrounding  the 
structure  was  varied  from  a  narrow  strip  only  10  ft.  wide  to  a  maximum  width  of 
90  ft.  Gne  group  of  structures  provided  buildi  ig  iength-to-width  ratios  from  1  to  co. 
The  bases  of  these  buildings  v  ere  20  by  20  ft.,  20  by  40  ft.,  and  20  ft.  by  co, 
calculations  being  made  for  each  for  three  wall  mass  thicknesses  (9,  36,  and  144  psf). 
A  second  group  consisted  of  two  square  buildings  20  by  20  ft.,  and  40  by  40  ft.,  again 
computed  for  wall  thickness  values  of  9,  36,  and  144  psf.  The  building  height  was 
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10  ft.  ond  the  detector  was  centrally  loccted  at  5  ft.  height  for  all  cases.  Computed 
values  of  the  direct  ond  scattered  radiation  for  the  two  groups  of  structures  are 
summcrized  in  Tables  5  and  6. 

In  certain  cases  the  values  of  direct  radiation  presented  in  Tables  5  and  6 
are  slightly  different  from  those  which  wouid  be  computed  using  the  data  of 
Figures  19  and  20  because  the  values  of  Gq  used  in  computing  these  results  varied 
slightly  from  the  refined  values  presented  in  Figure  20.  The  differences  in  direct 
radiation  caused  were  generally  less  than  10  per  cent  and  do  not  affect  the 
conclusions  presented  in  this  chapter. 

Shielding  calculation  of  the  Engineering  Manual  type  as  illustrated  in 
Figure  18,  for  this  series  of  single-story  structures  subjected  to  finite  rectangular 
fields  of  contamination  may  be  simplified  to  the  following  equation: 


Dose  =  direct  radiation  component  +  scattered  radiation  component 


D 

o 


B(X  ,h) 
e 


Gd(o.  ,h)  -  Gd(U*,h)'j  (1-SJ  ♦  Bos(Us,Xe)  [g  (»,  +G  <V> 


S  E 

W 


The  skyshine  component  is  omitted  from  this  equation  so  as  to  permit  direct 
comparison  with  the  NBS  Monte  Carlo  compilation  method.  The  NBS  method 
neglects  the  effect  of  the  atmosphere  and  assumes  that  all  radiation  enters  through 
the  detector  story  wails  of  the  structure  (i.e.,  there  are  no  radiation  components 
from  the  story  below  or  above).  The  direct  radiation  attenuation  factor,  B(Xe,5‘), 
was  taken  from  Figure  1 9b,  and  the  angular  distribution,  Gj(u  ,5')  and  G<j(u*,5'), 
generally  from  Figure  20b.  Two  complications,  however,  enter  into  the  calculation 
of  the  direct  radiation  components  as  the  building  width/iength  ratio  approaches 
zero.  The  first  arises  in  calculating  solid  angle  fractions  for  an  infinitely  long 
building  and  *he  second  is  introduced  by  the  tacitly  assumed  equivalence  of  circular 
and  rectangular  solid  angles  when  com,  ing  angular  distributions.  For  an  infinitely 
long  structure,  the  solid  angle  fraction  <j  equals  _  tan  “I  ( W  )r  where  W  is  the 

*  2Z 

structure  width  and  Z  is  the  detector  height.  As  for  the  second  complication,  one 
finds  that  in  differencing  G^O-^h)  for  a  series  of  buildings  with  fixed  width  and 
varying  length  and  a  fixed  finite  field  width,  the  value  of  A  Gd  first  decreases 
with  increasing  building  length  as  one  wouid  expect,  then  as  the  building  increases 
further  in  length,  A  Gd  increases.  As  building  length  increases  toward  infinity 
the  direct  contribution  and  hence  A  Gj  should  asymptotically  approcch  a  minimum 
at  a  length-tc-width  ratio  of  infinity  and  should  not  increase.  In  the  actual 
calculations  for  direct  radiation  for  a  given  source  field  and  building  width,  AGd 
was  plotted  versus  building  eccentricity  until  the  minimum  was  reached.  Values 
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beyond  this  minimum  were  extrapolated  with  the  extrapolated  curve  becoming 
tangent  to  o  constant  value  as  V-/L  -  0.  Thc^.  extrapolated  values  were  as  much 
as  10  per  cent  less  than  hie  calculated  —  value. 

The  NSb  Compilation  Prog^cm  for  finite  Fields  of  Contamination 

2q 

Eisenhas/er  has  recently  completed  a  program  to  compute  the  effect  of 
a  finite  field  of  contamination  using  transmission  and  scattering  coefficients  com¬ 
puted  by  the  Monte  Carla  method  by  Berger.^  The  advantages  of  this  type  of 
calculation  are,  first,  that  it  provide-  a  completely  different  method  of  calculation 
starting  with  different  assumptions  and  a  different  set  of  parameters  and  input  data, 
against  which  the  current  limited  stri>;  calculation  may  be  compared  and,  second, 
since  this  program  utilizes  transmission  coefficients  computed  by  the  Monte  Carlo 
method  it  may  truly  be  thought  of  as  a  number! cal  representation  of  the  cctual 
experiment.  This  code  has  been  designed  to  c  mpute  the  dose  contribution  to  a 
detector  located  behind  a  vertical  wall,  from  contamination  locc‘  *  in  a  rectangular 
area,  calk  1  o  source  pofh  (with  its  axis  parallel  and  perpendicular  to  the  wall) 
through  i  verficu'  rectangular  sec  den  .tie  wall.  Detector,  wail  section,  and 

source  patch  may  be  iocatea  at  any  arbitrary  position  in  relation  to  each  other, 
individual  cont.  but? on-,  from  each  source  patch  representing  differential  areas  of 
a  field  are  computed  a;,d  presented  separately.  Since  the  calculation  is  designed 
for  finite  fields  of  contamination,  air  scattering  and  attenuation  are  neglected. 

That  this  effect  is  negligible  for  finite  fields  of  the  size  of  interest  may  be  quickly 
demonstrated.  The  dcso  late  from  a  limited  circular  field  of  co^t amination  to  a 
detector  located  a  height  h  above  it  is: 


ft 

o 

D  *  2ttS  q  dp 

°  P 

h 


where  p 


toral  cross  section  ior  cir-  J/450  ft.  S.T.P. 

specific  irradiance  of  the  source  (raas/hr  curie)  at  1  ft. 

source  density  (curies/ sq.  ft.) 

detector  height  (it.) 

outer  siant  radius  of  the  field  (ft.) 
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vAieh  in  t is  absence  of  scattering  {all  collisions  are  assumed  absorption,  clearly 
placing  an  upper  limit  on  the  effect  of  the  atmosphere)  reduces  f? 


D  =  2a 


Soq 


— 

E 


E,G»H)- 


1 


(pR  )! 

°J 


where 


E 


1 


the  exponential  integral  of  the  first  kind 


if  the  ■c;««osoh*re  is  totally  neglected  p  is  zero  and  the  corresponding  equation  becomes: 


D  =  2.  Soqln  ^ 


The  maximum  error  introduced  by  neglecting  the  atmosphere  would  thus  become 


Error  = 


E,(ph)  -  £,(pRo) 
- 

In  £ 


Field  Radius  (ft.)  =  0 
Error  (%)  =  0 


10  20  30 

1.2  2.2  3.0 


50  70  100 

4.3  5.5  6.8 


The  dose  resulting  from  uncollided  photons  is  computed  in  Eisenhauer's 
program  by  first  determining  if  a  straight  line,  from  the  center  of  the  source  patch  ro 
the  detector,  passes  through  the  wall  area  (see  Figure  21).  If  it  does,  it  is  assumed 
that  all  the  contamination  in  the  rectangular  source  patch  is  concentrated  at  the 
center  of  that  ^tch.  The  dose  at  the  detector  is  computed  by  reducing  the  output  of 
that  contamination  by  the  distance  to  the  detector  squared  and  exponential  attenuation 
along  the  slant  distance  through  the  wall  (determined  by  the  straight  line  connecting 
the  source  patch  to  the  detector  location).  Scattered  radiation  is  computed  in  a 
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of  the  rec’mqirici  wall  crea  of  interest  and  then  drewing  a  second  line  from  this 
point  !o  the  defector  'location.  The  o.'gle  mode  by  these  two  lines  is  then 
determined  one  a  transmission  coefficient  for  scattered  radiation  is  interpolated 
from  the  TAonfe  CsHo  transmission  foci's  of  &erc?rr  The  radiation  emanating  from 
ihe  scarce  patch  is  then  reduced  by  the  radius  to  me  cer-ter  c?  th=  •*«!!  ore-  squared, 
the  transmission  vOeftTci cr-t  o°d  the  radius  to  the  detector  square". 


Clearly,  both  the  waii  oris  score-?  patch  grid  *?«»*  be  ct-:ose,?  with  cure  so 
that  furih-j.:  in  grid  size  will  yield  a  negligible  charge  in  resists, 

fxpipraicry  a.  eolations  indicate  that  for  a  detecsc-i  one  or  more  heights  behind 
the  wall,  locofed  at  the  wall  •mid-height,  a  square  wo'.*  grid  -,>ne-hc{*  I  height 
in  dimension  together  with  source  patches  of  a  similar  size  for  the  firsi  l  «c  wail 
h  .rights  av/cy  from  the  wcil  ana  twice  this  size  beyond,  aroduces  results  that  are 
within  i  per  cent  of  those  produced  b,  a  grid  one-fourth  as  large.  Eisenhauer^ 
estimates  that  this  program  produces  results  which  are  accurate  to  about  5  per  cent  if 
the  grid  sizes  are  properly  chosen. 


After  preliminary  evaluation  of  the  program  the  above-mentioned  grid  size 
was  selected  as  adequate  for  these  studies.  The  resultcnt  dose  contribution  from  o 
source  pctch  on  the  right  of  the  centerline  through  the  de'ector  perpendicular  to  the 
wall  are: 

1.  The  direct  contribution  —  the  dose  contribution  that  proceeds 
directly  to  the  de’.ector  without  scattering.  Since  the  detector 
is  located  at  th*.  vail  mid-height,  contribution  arises  through 
the  lower  half  of  the  wall  only  (direct). 

2.  7  ie  <vull-scatter  contribution  through  a  wall  slab  in  the  lower 
hnlf  of  the  wall  to  the  right  of  the  centerline  (wai!  scatter  lower). 

3.  The  .ymmetric  wall-scatter  contribution  through  a  wall  slab 
symmetrically  located  on  the  left  of  the  centerline  (symmetric 
scatter  low-.-'). 

4.  The  wall-scc'rer  contribution  through  a  wall  slab  in  the  upper 
half  of  the  wall  to  the  right  of  the  centerline  (wall  scatter  upper). 

5.  The  symmetric  wall  scatter  contribution  through  the  upper  wall 
slab  symmetrically  located  from  (4)  to  the  left  of  the  centerline 
(symmetric  scatter  upper). 

Although  the  balance  between  direct  end  scattered  radiation  is  set  by 
the  wall  thickness,  the  division  between  the  modes  of  scatter  penetration  seems 
quite  general  (at  least  for  the  cases  investigated,  between  one -quarter  and  four 
mean  free  path  thickness).  For  source  patcher  located  close  to  the  wall  the  wall- 
scatter  from  '.he  upper  half  of  the  wall  is  ger.eially  from  5  to  20  per  cent  of  wall 
scatter  from  the  lowc'  portion.  As  the  distance  from  the  source  patch  to  the  wall 
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increases,  this  r'ivisicn  becomes  more  neai'y  equcs,  «iih  values  generally  within 
10  per  .  snt  of  one  another  at  80  ft-  distance.  Scattered  radiation  arriving  through 
symmet..cclly  located  dctecto:s  generally  behoves  similarly  when  toe  source  patch 
end  wail  slob  ere  on  er  near  the  centerline.  As  either  the  slab  or  the  source  patch 
moves  off-center,  hc---.se ver.  ;;>>  symmetric  scatter  radiation  contribution  decreases 
mere  markedly.  Results  obteint  or  a  20  by  20  ft.  square  building  with  a  5  ft. 
detector  height  are  illustrated  in  figure  22  for  each  radiation  component.  The 
important  feature  ct  Figure  22  is  the  s'^ope  of  the  curves  cs  the  width  of  the 
contaminated  area  increases.  Ecch  component  of  transmitted  radiation  increases 
■*i  almost  the  same  way. 

The  NBS  Monte  Carlo  compilation  program  is  both  extremely  flexible 
and  easy  to  use  (though  summarizing  the  results  in  a  form  usable  for  this  study 
necessitated  much  manual  'abor  in  generating  the  required  number  of  "numerical 
experiments”  w:*h  which  the  cresent  method  c-f  calculation  could  be  evaluated. 

As  an  initial  stud-,  only  the  single  detector  height  or  5  ft.  was  chosen,  together 
with  c-  'taicfure  woi'i  height  of  10  ft.  Building  size  and  wall  thickness  were  then 
varied  and  the  data  zed  :w  t’-mt  the  finite  rectangular  limited  fir  kr 

of  contamination  surrounding  a  structure  could  be  cc'culo^d.  The  data  obtainc-d 
from  this  series  of  calculations  are  summarized  in  Table  7. 


The  two  methods  of  calculation  compared 

As  illustrated  earlier  in  tr>*s  chapter,  rhe  maximum  effect  of  air  scatter 
and  attenuation  on  the  radiation  originating  from  contaminated  areas  up  to  iOO  ft. 
wide  is  only  approximately  7  per  cent.  If  both  methods  of  calculation  produce 
the  same  results,  agreement  should  be  within  this  7  per  cent,  tbs-  engineering 
manual  results  being  lower  rr.cn  those  computed  by  rhe  NBS  Monte  Carlo  compilation 
method.  Figures  23,  24,  and  25  present  this  comparison  for  total  dose  a  function 
of  width  of  the  contaminated  field.  Figure  22  represent',  a  cose  in  which  the 
eccentricity  (W/L)  of  the  building  is  varied  from  zero  to  one  whh  all  other 
parameters  held  constant.  Figure  24  presents  a  series  of  caint.- lotions  in  w'-'icl-.  on.y 
the  wall  thickness  is  varied;  Figure  25  represents  two  cases  ider.fi col  except  for 
building  plan  area.  These  graphs  indicate  excellent  agreement  in  iota!  dose  os 
paedicted  by  the  two  methods.  The  mean  value  of  the  ratio  of  the  total  dose  as 
computed  by  the  NBS  methods  to  that  computed  by  the  Engineering  Manual  method 
I  for  all  cases  investigated  is  .98  (when  *he  ratio  is  computed  at  V'^h  =  2,3, 4,6,8, 

10,12,  14  for  each  case  investigated),  with  average  values  rancyng  from  .91  to 
1.04  for  any  one  structure  design,  i?  was  expected  that  this  ic'io  would  be  slightly 
greater  than  one  (because  atmospheric  effect:  were  neglected  in  the  NBS  method, 
rather  than  less  than  one.  This  minor  difference  mcbably  waf  caused  by  the  accuracy 
to  which  the  engineering  method  graphs  may  be  read.  Since  the  ratio  is  so  cl  se 
to  one  and  since  the  structure  investigated  was  varied  in  size,  wall  thickness,  and 


56 


CCNESCO  consultants  in  engineering  sC'e.xS 


—  co 

N  c, 
CM  o 


«Ti  o 

—  >c 

—  o 


—  *C 
MT  O- 

o  o 
o  o 


:  '  i? 

o  c 


8 


CO 

CC  fv 

I'M  O 
co  0 

—  m 

m 

CO 

o  <o 

o  o 

O'  ~wt 

'O 

CN 

o 

—  c 

o  o 

o  o 

3  * 

d> 

o  o 
co  co 
o  o 


CQ 

< 


o 

U 

£  o 
c  •- 


o  ~  ^ 

.=  O)  ^ 


i/i  6 

co  o 

Z  c 

o  8 


_D  ^ 

T3  <0 
P  Cl 


P 

_c 


o 

p 

-*— 

p 

TJ 


-O 

-a 

i 

U- 


c.  o 
E  3 
8  “ 
in  O 

°  o 

K  ai 

c  »- 

.2  o 

3  S 

-Q 


O  - — ' 

^  I 

v.  +- 
3  O 

4-  •  — 

U  TJ 
D  O 
i.  v- 

«*  > 
C  « 

o  E 

LT'- 

S.  • 


8  ^ 

o  «, 

4-  O 

o  o 
2  a. 

T3  C 
-a.2 

o  JB 


”D 

0) 


o 

u 

oo 


-C 

O) 

c 

c 

D 

o 


8 


o 

Mf- 


O 

co 


<3 


Q. 

c 

o 

o 


m 


-o 

o 

ct 


c 

p 

c 


E 

o 

u 


CM  O' 
CO  CM 
CM  O 


in  — 


~o 

p 

u  £ 
2  a 
O  £ 


CO 

8 

Mf  CO 

CO  CO 

CO  — 

CO  M3 

«o 

00 

O' 

o  o 

CO  mt 

O  CO 

rv 

CM 

o 

o 

o  o 

o  o 

o  o 

o 

o 

O  CO 

M  M3 

00  00 

co  rM 

fN 

M3  CO 

>c  «n 

D  CM 

00  M* 

o  o 

rv  cn 

•O  CO 

<2  CN 

Si  O 

o  o 
•  • 

O  O 

q  q 

q  q 

q  o 

m  cm 

^  m 

MJ-  — 

CN  O 

CN  cn 

GO  CD 

Tf 

o  o 

M3  CO 

Tj-  CN 

in  cn 

•  T 

o  o 

o  o 
•  • 

q  o 

o  o 
•  • 

c  o 
•  * 

O'  rx 

O  CO 

co  co 

—  N- 

CO  «“! 

tM  co 

O  M3 

o  — 

in  co 

o  o 

T*  CN 

CO  CM 

N-  - — 

-  o 

•  • 

O  o 
•  • 

o  o 
•  • 

q  (?- 

O  O 
•  • 

q  o 

m  rx 

4-  cr> 

rv  rs. 

—  co 

O'  — 

-  O' 

CM  CO 

—  - — 

if  CM 

o  o 

CO  CM 

C*i  *— 

00  - — 

-  o 

o  o 
•  • 

o  o 
•  • 

o  o 

q  q 

o  o 

—  --O 

t  ™ 

m  o 

•—  CN 

O'  m 

co  -e- 

oo  o» 

CO  CM 

o  o 

CM 

CN  — 

CN  O 

3  O 

l  • 

o  o 
•  • 

o  o 

q  q 

o  o 

q  o 

M3  CO 

0  M3 

o  o 

o  o 

M3  CO 

00  !n 

CO  Mf 

*  o 

CN  .*— 

o  o 

—  O 

—  o 

—  o 

D  O 
•  • 

O  10 

o  o 

a  o 
•  • 

o  o 

o  q 

fl) 

ul  o 

5  £ 


"O 

P 

U  £ 

d)  *■ 

£:  o 

5  £ 


•o 

a 

L. 

",  P 


"O 

p 


o 

U 

1.0 


o 

u 

in 


"O 

p 


o 

u 

1/1 


p 

—  c 

"5  -X  t*. 

.2  H 


o-  o- 


MO  MO 
CO  CO 


33 


MO  MO 
CO  CO 


MO  MO 
00  CO 


Ml  MO 
00  (O 


O) 

c 


3 

CO 


a 

O 

o 

o 

o 

CN 

CN 

CM 

N" 

8 

X 

,-< 

X 

X 

X 

X 

o 

O 

o 

o 

o 

o 

CN 

CN 

CM 

CN 

CN 

57 


/ 


Normalised  so  as  to  make  the  response  3  ft.  above  or  infinite  field  of  contominati 
1.0  if  *he  effects  of  the  atmosphere  were  included. 


0  10  20  30  40  50  60  70 

Field  Width  Wc(ft) 


TOTAL  GROUND  DOSE  FOR  BUILDINGS  OF  DIFFERENT  SIZE 
(36  p$f  Exterior  Wall  10  Feet  High,  Detector  Height  5  Feet) 

Figure  25 
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eccentricity,  it  may  be  safely  concluded  therefore  that  at  first-floor  locations  as 
typified  by  detector  height  of  5  ft.,  the  current  method  of  computing  the  effect 
of  finite  rectangular  fields  of  contamination  is  adequate. 

Unfortunately  there  was  not  enough  time  to  generate  a  similar  group  of 
Jata  for  comparison  at  other  detector  heights  (for  example,  in  bcsement  areas  and 
upper-floor  locations)  in  order  to  make  the  conclusions  more  general.  In  an  attempt 
to  estimate  the  effects  at  other  detector  locations,  the  scatter  and  direct  dose 
compared  component  by  component.  A  surprising  result  was  uncovered:  although 
the  ratio  of  either  direct  or  scattered  dose  as  computed  by  the  NB5  Monte  Carlo 
Method  to  that  computed  by  the  engineering  method  was  very  nearly  constant  as 
a  function  of  field  size  (significant  variations  from  the  mean  value  occurred  only 
within  the  first  15  ft.  of  field  size),  this  ratio  varied  significantly  from  one. 

Table  8  presents  a  summary  of  these  data  computed,  as  before,  for  each  case  investi¬ 
gated  at  field  widths  of  5,  10,  15,  20,  30,  40,  50,  and  70  ft. 


TABLE  8  i 

Mean  ratios  of  scatter  and  direct  radiation  as  computed 
by  the  engineering  method  (10  ft.  structure  height, 

5  ft.  detector  height) 


5  ! 


_ 

rati 

os 

W 

Parameter 

varied 

Structure 

(ft.) 

Wall  thickness 

(psf) 

Direct 

Scat. 

Old 

Indicated 

20  x  20 

9 

1.13 

.513 

.178 

.11 

Wall  thickness 

20  x  20 

36 

1.45 

.751 

.475 

.34 

20  x  20 

144 

1.92 

.846 

.828 

.70 

20  x  20 

36 

1.45 

.751 

.475 

.34 

Eccentricity 

20  x  40 

36 

1.39 

.714 

.475 

.33 

20  x  oo 

36 

1.08 

.800 

.475 

.37 

20  x  20 

36 

1.45 

.751 

.475 

.34 

Size 

40  x  40 

36 

1.48 

.541 

.475 

.27 

Sizable  variations  between  the  ratio  of  direct  and  scattered  radiation  occur 
as  a  function  of  wall  thickness,  eccentricity,  and  structure  size.  Because  data 
available  are  limited,  it  is  difficult  either  to  pin  down  the  reasons  for  this  variation 
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or  to  recommend  a  change.  Since  the  ratios  of  total  radiation  computed  by  both 
methods  were  so  nearly  one,  however,  the  ratio  of  scattered  to  total  radiation, 
Sw(Xe),  used  in  the  current  methodology  becomes  suspect.  The  values  of  this 
parameter  required  to  make  the  ratios  of  Table  8  equal  to  one  are  presented  in  that 
table  and  in  Figure  26,  together  with  the  values  currently  used.  Because  of  this 
discrepancy  we  must  therefore  caution  that  at  this  time  the  conclusion  that  the 
engineering  method  of  computation  is  valid  for  limited  fields  of  contamination  may 
be  made  only  for  locations  on  the  first  floor  of  a  structure.  Other  locations,  where 
the  relative  importance  of  direct  or  scattered  radiation  is  considerably  different 
(as  in  basement  areas  or  upper-floor  locations  where  the  field  is  shadowed  by  a 
thick  floor)  must  be  investigated  further. 
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Wall  Thickneu  X#  (p$f) 


FRACTION  OF  TOTAL  RADIATION  THAT  IS  SCATTERED 
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CHAPTER  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  study  we  have  examined  the  new  method  ^proposed  for  computing 
the  attenuation  of  "in  and  down"  scattered  radiation  and  the  current  method  of 
computing  the  effects  of  finite  fields  of  contamination  so  that  we  might  recommend 
improvements  and  up-dating  for  OCD  documents  OCD-TR-20,  Volumes  1  and  2, 
"Shelter  Design  and  Analysis." 


in  and  Down  Scattered  Radiation 


The  new  barrier  factor  B'0(X/  u  )  computed  for  in  and  down  scattered  radiation, 
based  on  the  assumption  of  Compton  first  collision  scatter  in  the  structure  walls,  has 
been  found  to  provide  significantly  better  agreement  with  existing  experimental 
measurements.  Comparing  this  factor  with  approximately  350  experimental  measure¬ 
ments,  we  found  that  the  new  method  was  generally  conservative,  predicting  dose 
rates  from  ground-based  sources  of  contamination  18.5  per  cent  higher  than  those 
measured,  with  a  standard  deviation  of  13.5  per  cent. 

It  is  recommended  that  this  new  barrier  factor  B'o  (X,u  ),  reduced  by 
18-1/2  per  cent,  replace  the  currently  utilized  Bq(X)  to  improve  the  accuracy  of 
the  calculation  method.  In  Chapter  3  of  this  report  we  provide  detailed  recommendations 
on  the  steps  required  to  effect  this  change  in  both  the  Engineering  Manual, 

OCD  TR-20,  Vol.  1,  and  the  equivalent  building  method,  OCD-TR-20,  Voi.  2. 

Finite  Fields  of  Contamination 

A  preliminary  investigation  to  determine  the  effects  of  finite  fields  of 
contamination  for  detectors  on  "first-floor"  locations  indicates  excellent  agreement 
between  total  doses  when  calculated  by  the  two  different  methods,  'vloior  variations 
in  the  scatter  to  direct  ratio  were  obtained  by  the  two  methods.  There  is,  therefore, 
serious  doubt  as  to  the  accuracy  of  the  recommended  procedures  for  positions  where 
the  scatter  to  direct  ratio  is  significantly  different  from  that  investigated  (e.g.,  in 
basement  or  upper-story  locations  where  the  detector  is  shadowed  by  a  heavy  floor 
slab).  It  is  recommended  that  the  effects  of  finite  fields  of  contamination  in  areas 
where  direct  radiation  does  not  provide  a  major  dose  component  be  studied  further 
in  order  to  develop  a  simplified  method  of  calculation  that  will  be  accurate  in  all 
locations. 
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General 


The  current  program  has  reviewed  most  of  the  experimental  and 
analytical  efforts  devoted  to  shielaVtg  from  contamination  in  a  "fallout  geometry." 
It  has  become  evident  during  this  review  that,  in  addition  to  the  two  areas  of 
inquiry  that  are  covered  in  this  report,  several  additional  areas  should  be  investi¬ 
gated.  These  are:  the  ratio  of  scattered  to  direct  radiation,  particularly  as 
applied  to  finite  fields  of  contamination  in  below-ground  and  upper-story  shelter 
locations;  a  preliminary  investigation  of  simplified  methods  for  computing  the 
effects  of  inhomogeneous  floor  and  wall  barriers;  and  a  preliminary  study  of  the 
effects  of  roof  wash  down  (or  other  decontamination  methods)  and  how  these 
effects  may  be  simply  expressed  (perhaps  as  an  additional  equivalent  roof  thickness) 
in  each  of  the  present  methods  of  caltulation.  The  latter  program  is  suggested  so 
that  this  procedure  (of  decontamination)  may  become  useful  in  analyzing  new 
building  construction  when  other  slanting  or  design  techniques  cannot  be  used. 
Additional  effort  is  required  in  each  of  these  areas  to  ensure  that  the  appropriate 
methods  of  analysis  currently  used  or  proposed  are  accurate,  simple,  and  easily 
usable  so  that  they  may  be  included  with  confidence  in  OCD  publications. 
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